Replicating financial market dynamics with a simple 
self-organized critical lattice model. 



B. Dupoyet a ' b , H.R. Fiebig a ' c , D.P. Musgrove a,c 

"•Florida International University, Miami, Florida 33199, USA 
b Department of Finance 
c Department of Physics 



Abstract 

We explore a simple lattice field model intended to describe statistical properties of 
high frequency financial markets. The model is relevant in the cross-disciplinary area of 
econophysics. Its signature feature is the emergence of a self-organized critical state. This 
implies scale invariance of the model, without tuning parameters. Prominent results of 
our simulation are time series of gains, prices, volatility, and gains frequency distributions, 
which all compare favorably to features of historical market data. Applying a standard 
GARCH(1,1) fit to the lattice model gives results that are almost indistinguishable from 
historical NASDAQ data. 

Keywords: Econophysics, Financial markets, Statistical field theory, Self-organized 
criticality 
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1. Introduction 

From a reductionist perspective the statistical physics of a large number of dynamical 
systems in nature originates from nonlinear processes at the microscopic level. In many 
cases this leads to phenomena characterized in the literature by the terms of chaos, 
complexity, fractal geometry, and criticality. These scenarios are quite ubiquitous, thus 
not limited to basic physical systems, where turbulence comes to mind for example, but 
also to applications in biology, geology, social networks, economic systems, and finance, 
to name a few The literature on the subject is prodigious. 

Our current interest in the subject stems from a recent simulation of financial market 
dynamics At the root of that study is a microscopic model based on the principal 
of gauge invariance, assuming that one of the key mechanisms of trader behavior is 
independent of any scale (currency unit, for example) used in the market transactions 
In technical terms the model is a quantum field theory based on the gauge group 
G = R + , the dilation group, which implies scale invariance of the market model with 
respect to ordinary multiplication of prices with positive real numbers. The quantum 
aspect of that model implements the empirical observation that arbitrage opportunities, 
i.e. realizing a profit via transactions in different markets, vanish quickly because of 
market dynamics. 
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At this stage, the model does not provide a mechanism for describing a complex 
system, as it should, given the empirical evidence. The distribution of market returns, if 
analyzed appropriately Q , exhibits fat tails (probabilities larger than Gaussian) the likes 
of which are observed in many high frequency financial markets. However, this is not a 
consequence of the intrinsic dynamics of the model. In order to remedy this situation, 
in the present article, we study an abridged model with local interactions that lead to 
a self-organized complex market model. Although our goal is to eventually combine the 
gauge model with features of the abridged model discussed here, the latter, despite its 
simplicity, produces salient characteristics of actual financial markets surprisingly well. 
Among those are the semblance of return time series, returns frequency distributions, 
and the nature of volatility. The market volatility, in particular, is a subject of intense 
research 0, [1, H . 

These features are promising enough to study this simple model in its own right, 
which is the subject of this work. Along a one-dimensional lattice representing discrete 
time, the field on the sites are interpreted as returns in a model market. An updating 
algorithm is then applied which is loosely fashioned after the well-known proposal by Bak, 
Tang and Wiesenfeld (BTW) Q, see also @, i, ED, Gil ■ The idea is to develop a market 
model that is driven by microscopic entities, say traders, such that their interaction leaves 
the lattice field in a self-organized critical (SOC) state. In a critical state, among other 
things, long-range correlations of suitable observables lead to power-law behavior with 
respect to scaling transformations. Self organization means that the system is driven 
to criticality without fine-tuning any external parameters, i.e. solely by its intrinsic 
dynamics. 

It is generally realized that financial markets, being prime examples of social systems, 



exhibit SOC [11|, |12| . Il3l . Il4 Il5| . However, to the best of our knowledge, attempts to 



model those from a microscopic point of view are rare 16| . In 17|, [l8[ percolation clusters 
act as investors. Assigning random percolation probabilities, power law behavior is found 
for the usual observables derived from stock market prices [19| . Another ex amp le of such 
an attempt close to the BTW evolution model is the work of Bartolozzi et al [20( . Though 
close to our work in spirit, significant differences exist in terms of the updating strategy 
and interpretation. Our implementation of our lattice field model will produce price time 
series, returns and their distributions, volatility time series and their clustering features 
that are all strikingly similar to historical market data. 

Why would one like to have a market model in the first place? After all there are 
myriads of historical data being collected every day. A good reason is the fact that all of 
the collected data are merely instances of a random draw from some probability distri- 
bution, just like one throw of a dice gives only one number, hiding the statistics behind 
it. A stochastic model on the other hand will enable us to study any number of market 
instances, and collect ensembles in the language of statistical physics. Observables, as 
averages endowed with errors, could be computed. Ultimately, a successful model could 
provide probability distributions for future prices, and thus be an invaluable tool for risk 
analysis, and the like. 



2. Lattice model 



We consider the simplest lattice market model conceivable, a one-dimensional chain of 
n+ 1 sites with labels j = . . . n, where j indicates discrete time t = jA in steps of some 
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arbitrary unit A. The sites are populated with a real- valued field r with components 
rj G K. When compared to [2( , we use a minimalistic topology, and also ignore the gauge 
field living on the lattice links thus forgoing the minimization of arbitrage. As it turned 
out it is essential to interpret the field components rj directly as investment returns. The 
returns are defined as 

r s = log^/^.r) (1) 

where $j = Pj/C is the price of an investment instrument, such as a stock or index fund 
for example, and C is a unit (currency, shares, etc). The continuum version of ([1]) can 
be surmised from taking the limit A — ¥ in 

k>g($(f)/$(t - A)) = A^ log *(t) + 0(A 2 ) , (2) 

where <&(i) = at t = jA. 

In order to endow the field r with dynamics we find inspiration in the popular evo- 
lutionary model by Bak and Sneppen 0, [l], ■ In that context, the field components 
are fitness values, say fj € [0, 1], assigned to the sites of a lattice. The updating process 
consists in finding the site j s with fj s = minj/j : j = 0...n}, i.e. the least adapted 
species. Then fj s and the values fj s ±i of the two next neighbors are replaced with uni- 
formly distributed random numbers from [0,1]. This prescription, when iterated many 
times s = 0, 1 ... oo, leads to a stationary state of the lattice field where a single pertur- 
bation can lead to a burst of activity, called an avalanche. The frequency distribution of 
avalanche sizes is found to follow a power law. A power law is a signature feature of a 
critical state. Since no tuning of a model parameter is needed, the phenomenon is known 
as self-organized criticality (SOC). The model is very robust in the sense that changing 
the updating prescription, within reasonable bounds, will still lead to SOC. A rigorous 
discussion, containing analytical results, may be found in fioj . 

In the context of the financial market model we adopt a modified version of Bak's 
updating prescription. We select periodic boundary conditions with period n + 1, such 
that r n+ \ = ro and r_i = r n . In terms of the returns rj, we define 

Vj = rjirj+i-r^i) (3) 

Vj = N (4) 

and call 

V = max{Vj ■: j = 0...n} (5) 

the signal of the field configuration r. The updating strategy then proceeds with finding 
a site j s fronfj] 

j,e{j = 0...n:V s =V}, (6) 
and then replacing the returns on sites j s and its two neighbors according to 

rj s xq and r,- s ±i 4- x±\ (7) 

where xo and x±i are three random numbers drawn from a normal distribution p(x) oc 
cxp(— x 2 /2w) while enforcing the constraint X-i + xq + x + \ — 0. The variance w is a 
parameter. Figure [1] illustrates the situation. 



1 There is at least one, and almost always only one element in this set. 
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Figure 1: Illustration of the geometry of the lattice model and the label scheme for the sites. Periodic 
boundary conditions rj+n+i = rj are implemented. Updating is done on the field component with signal 
V, and its two next neighbors. 



Within reason, we have experimented with numerous alternative definitions for the 
signal, defined through (|3I4I5[) . Overall, it appears that Bak updating is very robust and 
SOC is easily achieved. However, the choice (|3I4I5|) proved to best match the stylized 
features of historical financial market data. 

Although (|3I4I5|) were mostly selected on empirical grounds, in retrospect, more mo- 
tivation may be provided: Note that ([3]) is just a discretized version of 

■<*> = r(i)2A ^ir - a s r "> 2 < 8 > 

where r(t), in view of ([T]) and ©, has been identified with 

r(t) = A^ log *(t). (9) 
In finance, an established approach is to treat the returns r(t) as a stochastic process 



211 ] . Returns are typically modeled by a (generalized) Wiener process, i.e. assuming 



normal distributed random variables with a time dependent (random) variance 

W(t)=E[r(tf}-E[r(t)} 2 , (10) 

where £[•••] indicates the stochastic expected value. Discretized versions of the time 
derivatives of the two terms are 

jE[r{tf]=2E[r{t)r'{t)] ~ < rj (r i+1 - r^A" 1 (11) 
± t E[r{t)Y = 2E[r{t)]E[r'{t)] ~ ( rj )<r i+1 - r^A" 1 . (12) 

On the right-hand sides, we have changed the notation for the expectation value from 
E[- ■ ■} to (•••), the latter indicating the averages for lattice generated returns. In the 
current simulation, standard stochastic financial modeling dictates that (rj) be indepen- 
dent of time. This implies that the discretized part of equation (12) is exactly equal to 
zero. 

Therefore, the dynamics of the lattice model is driven by eliminating extreme, sudden, 
changes of the variance on the returns deemed 'unfit' in the spirit of [9(. Combining (|10[) 
and (|llll2p the lattice version of the latter turns out to be 

dWit) , , , 

— — ~ [cov(rj,r i+ i) - cov(rj,rj_i)]A , (13) 
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where cov(r a ,rp) = (r a rp) — {r a )(rp) is the covariance of the two random variables. 
Hence, extreme changes of the covariance of returns between adjacent time slices are 
discouraged as part of the dynamics of the market model. 

As a side remark we comment on the use of the two-step time derivative to approx- 
imate r'(t) in (TTTj) and (|T2"jl as opposed to employing one-step forward (rj+i — r^A -1 , 
or backward (rj — rj_i)A , discretizations. In both cases the interaction ([3]), which 
drives the lattice dynamics, would mutate from a proper next-nearest-neighbor cou- 
pling to a term dominated by self- interactions cx r?. This does not lead to a sensi- 
ble physical system, and demonstrably gives absurd results in a numerical simulation. 
The corresponding approximations to the time derivative of the variance dW/dt are 
~ [cov(rj + i, rj) — var(rj)]2A _1 and ~ [var(rj) — cov(rj, r :) _i)]2A~ 1 , respectively. Those 
are variance driven and, equally, not suitable to define market dynamics. 

Finally, we should caution that the purpose of the above exposition is only to provide 
some motivation. In fact we should expect that the returns time series has fractal geom- 
etry. This makes perfect sense on a discretized time lattice, with the usual implications 



22t l23j , but it certainly discourages using the concept of continuous time derivatives on 



everyday manifestations of market data. In this vein the lattice model, in combination 
with next-nearest-neighbor interactions, appears to be a rational approach. 



3. Simulation 

For good measure, we have chosen a lattice with n — 780 time steps. Starting with a 
random lattice field r we show in Fig.[5]an example of the updating evolution of the signal 
V = max{|rj(rj + i — : j = . . .n} versus the simulation 'time' s. Clearly visible is 

a significant drop of the lower envelope of V(s) as s approaches ~ 10000. Beyond that 
the distribution of signal dots appears to have stabilized at around « 20000, at least this 
is the visual impression. 

Following [lOj we take a closer look at the envelope and define the 'gap' function 

G(x) = min{V(s) : s e N U {0} and s < x} with x € K + U {0} , (14) 

which is meant to trace the lower envelope of the signal. By construction, G{x) is 
a decreasing piecewise constant function with discontinuities at certain discrete values 
Xk , k G N. Adding xq — and assuming an ordered sequence the length of a plateau is 
A/c = Xk — Xk-i, and its height is G(xk-i). By definition, we say that an avalanche of 
length Afc starts at Xk-i and ends at Xk- At s — Xk-i all lattice sites have local signals 
Vj < G(xk-%), see ([5]). As long as the avalanche lasts, there is at least one lattice site 
with a local signal larger than G(xk-i) and thus the updating activity continues until 
s = Xk- Since the gap function is decreasing and bounded from below by zero it will 
eventually approach a constant lim ;E _ s . 00 G(x) — Gc- In this regime the avalanche size 
diverges and the system has reached the desired state of criticality [lj| • An example 
of G(x) from our simulation is shown in Fig. [3] It corresponds to the data of Fig. [5] 
but up to much larger simulation times. The evidence points to a critical value Gc with 
< Gc iS 0.01. Note that Gc = is a possibility, we do not know if it is realized. Unlike 
for the evolution model [l(| no analytic results are available. Strictly speaking, the above 
narrative applies to the thermodynamic limit, implying a lattice with infinitely many 
sites. Thus, in principle, we expect finite-size effects to afflict our simulation. However, 
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Figure 2: Updating evolution of the signal V versus the simulation time s. The signal at s = is due to 
a random start of the lattice field configuration. (Its value is = 4.74.) 

because subsequent results are very sensible we don't expect those to be an obstacle to 
practical application of this model. 

A signature feature of a critical system is scale invariance, implying power law behav- 
ior of certain quantities. Again following flpj . we display in Fig. 0] the frequency distri- 
bution of the avalanche sizes AN/ AA where AA is a binning interval for the avalanche 
sizes and AA is the count of avalanches within that interval. We have used 10000 bins 
with a binning interval of A A = 1. The data points come from an ensemble average 
over 2000 independent lattice simulations with 2 x 10 6 update steps each. This allows 
one to calculate statistical errors, also displayed in Fig. [U A power law behavior is be- 
yond doubt. A least-% 2 fit including data points in the interval 10 1 < A < 10 2 gives 
AA/AA = 301 A -139 . Integrating this gives a total number of A « 780 avalanches out 
of a run with 2 x 10 6 updates. However, the integrated average A is clearly less than the 
actual avalanche count per simulation due to the fact that the averaging operation gives 
rise to counts less than unity, which is realistically unfeasible. An analysis over many 
simulations gives a more accurate count average where 780 < A < 1000 for a single run 
with 2-million updates. Hence long avalanches clearly are key to the updating process. 

In Fig. [5]we display the activity of the lattice sites during an updating sequence. By 
definition, H is the number of times a lattice site j has been visited by an updating 
hit since the start of the simulation. Initially, with random Tj assigned to the sites the 
signal V tends to visit each site with comparable probability, leading to a flat activity 
plot. The plot in Fig. [5] reflects the activity after 20000 updates. A glance at Fig. [2 and 
Fig. O reveals that this corresponds to the onset of criticality. The peak-like structures 
in the activity plot give evidence of the emergence of avalanches with larger sizes. 

A closely related plot shown in Fig. [5] shows a zoom window on the simulation time 
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Figure 5: Plot of the activity across the lattice sites after 20000 update steps. The appearance of peaks 
signals the emergence of avalanche dynamics. 
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Figure 6: A zoom window on the updating evolution. Each plot symbol (-) at site j and simulation time 
s indicates an updating hit. 
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dependence of the updated sites around j ~ 420 and s ps 10 6 . Each plot symbol in 
Fig. [6] indicates an updating event of site j at simulation time s. Disconnected sets of 
dots clearly show the presence of avalanches. By visual inspection the curve is fractal in 
nature, a signature feature of complexity. In the context of the evolution model [s] the 
term 'punctuated equilibrium' has been used to describe a similar observation. 

Finally, we have used an entropy-like quantity to monitor the approach of the field 
towards a complex state. Using the exponentiated returns Rj = exp(rj) define 

1 " 

Then Fig. [7j in which is displayed S versus the simulation time s, shows that the initially 
random system becomes organized, again at around s = 20000. From then on the 
information content of the field has stabilized, as indicated by bounded fluctuations of S 
between 10 _1 and 10~ 2 . 
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Figure 7: Entropy S of the exponentiated returns Rj = exp(rj) as a function of the simulation time s. 



4. Results 

As mentioned above, with our lattice geometry and size, it takes at least 20000 
updating hits to reach a critical state. The lattice has n = 780 time intervals, and we 
use w = 1. To obtain the results discussed in this section we have used 4 x 10 6 initial 
updates before collecting field configurations from independent simulations. 

4-1- Gains distribution 

The advantage of working with a stochastic model is that observables can be estimated 
from ensembles, i.e. multiple realizations of market time series. Thus we can take a look 
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Figure 8: Lattice returns (gains) distribution. Ac/Ar is the number of returns of size r at a binning 
interval of Ar = 0.05. 
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Figure 9: Historical NASDAQ returns (gains) distribution (open circles, H) compared to the results of 
the lattice simulation (filled circles, L). The dashed lines arc Gaussian fits to the center and the tails, 
respectively, of the lattice data. 
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at the gains distribution, defined as the frequency plot of returns against the size of the 
returns. The lattice data displayed in Fig. [8] are averages from 10000 simulations, taking 
one sample after 4 x 10 6 updates each. The binning interval for the returns is Ar = 0.05 
and the number of counts per interval Ac/Ar is normalized such that the total number 
of counts is n = 780. Statistical errors are visible at the tails of the distribution. A 
well-known feature of gains distributions is that they exhibit 'fat tails', meaning that for 
extreme values of r the distribution is considerably enhanced over a normal (Gaussian) 
distribution. In [2|, where we studied the effects of arbitrage using a local gauge field, 
this feature was only obtained after analyzing the returns performing a weighted time 
average of past returns. In the present model the fat tails distribution emerges naturally. 
Since the returns field develops into a critical state the resulting scale invariance implies 
loss of memory of past holding patterns. 

The same gains distribution is displayed in Fig. |H] along with historical market data 
from the NASDAQ index compiled from minute data between 2005-Aug-26 and 2008- 



Aug-25 (2J. Scale factors as in AX = 2.4 x 10~ 3 Ar and AN/ AX = 1.1 x 10 5 Ac/Ar 



have been applied to match the historical data. The dashed lines correspond to Gaussian 
distributions fit to the lattice data in the center (7 points) and the tails (2 x 38 points), 
respectively. The Gaussian distribution is clearly ruled out. We observe that the lattice 
model accounts remarkably well to the empirical gains distribution over many orders of 
magnitude. 



4-2. Time series 

In Figs. [TU1[T3l we show four sets of sample time series from the lattice and selected 
historical data 24] from the NASDAQ index. The latter are included to demonstrate 
that the lattice model has the 'distinctive air' of a real market. It goes without saying 
that a statistical model is at a loss predicting time series, nor can one expect that a single 
model will describe subtle details of stochastic features of every market. What we are 
looking for has been very eloquently layed out in [22| by Mandelbrot. In the Introduction 
we read: "It is worth noting that fully fleshed-out and detailed pictures ... put a heavy 
premium on the ability of the eye to recognize patterns that existing analytic techniques 
were not designed to identify or assess." It is in this spirit that Figs. IT0l - fT3l are presented. 

The NASDAQ data from [24{ are at minute intervals, which naively translates to A = 
60s for the parameter introduced at the beginning of Sect. [2J There are discontinuities 
by end-of-the-day and over-the-weekend interruptions in the time series. Ignoring this, 
we have randomly picked time series of length n from the historical data. However, the 
situation is more complicated. Assuming that the time series has fractal geometry, and 
thus is devoid of a scale, the choice of any A would be equally valid. We are not in 
a position to pursue this issue here, but will rather live with the above naive choice, 
adopting the scaling argument. 

As for Figs. HUHT^l the three panels on the left column display, top to bottom, for 
j = . . . n the returns Tj , the prices pj and the volatility Vj of returns derived from the 
lattice model. While the returns come directly from the simulation, the price time series 
is obtained by integrating ^ 

1 I* 

$(i) = *(*„) exp[- / dt'r(t')]. (16) 
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The discrete version of (fl"6)) is equivalent to rewriting ([TJ as a recursion 



Pj =Pi-iexp(r i ) 



(17) 



with pj — QjC and initial condition pq. Finally, the time series of the variance, or 
volatility in financial terms, is computed from just three time slices through 



Note that A,po and C are all trivial parameters (units, etc) in the sense that they 
have no effect on the simulation. The same is true for the variance w used to draw 
random numbers, from p(x) oc exp(— x 2 /2w), while updating the lattice, see ([7]). Since 
the signal ([5]) of the held configuration only involves a comparison (max) of returns the 
effect of changing w will be a rescaling of the field. Specifically 



Therefore, doing an entire simulation at only one fixed, arbitrary, choice for w is suf- 
ficient. In this sense the model has no adjustable parameters. Observables then scale 
accordingly, for example Vj — » Xvj . The price time series (|17p behaves in a less straightfor- 
ward, though well defined, manner because the scaling factor appears with the argument 
of an exponential function. The plots in Figs. [T0HT31 were consistently produced with 
w = 1, A = 2 x 1CP 5 , and po = 1950. The resulting scales are similar to those of the 
selected historical NASDAQ data. Aside from lateral shifts of the price viewing windows, 
comparable scales in all panels of Figs. [TUHT51 are identical. 

The four data sets in Figs. ITUHT31 were chosen somewhat randomly, except that they 
were paired to point out common features appearing in both the lattice and historical 
sets. For example, in Fig. [10] there is a clearly visible time interval around j « 330 
where the historical returns data exhibits a region of quenched volatility. In financial 
market data the phenomenon of volatility clustering is well known and even exploited as 
an industry standard modeling technique [21] . Periods of small/large volatility are often 
separated by price shocks that manifest themselves in large spikes of the volatility. This 
is clearly a common feature of both the historical and the lattice data sets of Fig. [TOl 
best visible in the v panels. The price evolutions are similar, but this is accidental. The 
lattice model generates up or down markets with equal probability. (Although this can 
easily be changed by modifying the updating algorithm.) 

The remaining sets exhibit the same characteristics. In Set 2, Fig. [TTJ we have paired 
lattice and historical data that share a small overall volatility. This demonstrates that the 
lattice model is capable of producing quiescent market periods as well. The price times 
series look strikingly similar, which is fallacious, the statistical model is not capable 
of making predictions. The subsequent Sets 3 and 4 provide additional evidence that 
the lattice model is able to emulate stylized featured of real markets. Again, volatility 
clusters are clearly present separated by spikes of various sizes. Similarly, Set 4 exhibits 
fairly active volatility patterns for both the lattice and the historical data. 

4-. 3. Financial market dynamics 

We now turn to the ability of our model to capture some of the well-known dynamics 
observed in financial markets. One of the most important characteristics of a financial 




(18) 



w — > Xw then rj — > V~\rj . 



(19) 
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Figure 10: Set 1: Lattice generated times series of returns r, prices p, and volatility v are shown in the 
left column. The three panels on the right show selected historical r, p, v times series from the NASDAQ 
index. 



Set 1 L 



Set 1 N 



a 3.23E-08(1.01E-08) [3.214841] l.HE-08(5.35E-09) [2.079218] 
ai 0.043332(0.008898)[4.870008] 0.130702(0.019696) [6.636112] 
Pi 0.891148(0.028338) [31.44721] 0.859163(0.026981) [31.84381] 



Table 1: GARCH fit parameters of Set 1, FIG. 1101 for the lattice L and historical N data. 
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Figure 11: Set 2: See caption of Fig, 1101 



Set 2 L 



Set 2 N 



a 4.14E-09(1.61E-09) [2.564338] 4.20E-09(1.07E-09) [3.916649] 
ai 0.022450(0.005070)[4.428308] 0.031982(0.007263) [4.403283] 
ft 0.966177(0.007889) [122.4713] 0.951740(0.010025) [94.93278] 



Table 2: GARCH fit parameters of Set 2, FIG.[TT] for the lattice L and historical N data. 
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Figure 12: Set 3: See caption of Fig. \W\ 



Set 3 L 



Set 3 N 



a 6.38E-09(1.46E-09) [4.360360] 5. 87E-09(2.02E-09) [2.907091] 
ai 0.015515(0.003301)[4.700237] 0.022312(0.008688)[2. 568121] 
ft 0.972397(0.004485)[216.8057] 0.956839(0.014609) [65.49589] 



Table 3: GARCH fit parameters of Set 3, FIG.Q2] for the lattice L and historical N data. 
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Figure 13: Set 4: See caption of Fig. \W\ 



Set 4 L 



Set 4 N 



a 2. 76E-08(5.42E-09) [5.092728] 3.73E-09(2.70E-09) [1.378933] 
ai 0.045674(0.007692) [5.937874] 0.069064(0.005237) [13.18884] 
ft 0.911153(0.013877) [65.65998] 0.936347(0.003615) [259.0275] 



Table 4: GARCH fit parameters of Set 4, FIG.Q2] for the lattice L and historical N data. 
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time series is its volatility, and more importantly, how the volatility evolves over time. 
Most financial time series exhibit time- varying volatility clustering, which means that 
periods of large swings tend to be followed by periods of large swings, while periods 
of calm tend to be followed by periods of calm. These dynamics can be modeled by 



the Auto Regressive Conditional Heteroskedasticity (ARCH) model of Engle [21( and by 
its generalized version, the Generalized Auto Regressive Conditional Heteroskedasticity 
(GARCH) model of Bollerslev (H- 

The fact that these specifications imply that a large shock on average tends to be 
followed by another large shock means that the resulting distribution of returns will 
exhibit 'fat tails' or higher-than- Gaussian probability masses in the extreme regions. 
Similarly, the fact that these specifications also imply that a small shock on average tends 
to be followed by another small shock means that the resulting distribution of returns will 
exhibit a higher-than-Gaussian probability mass around the origin, see Fig. |H1 ARCH 
and GARCH models have grown into incredibly popular tools as they are able to replicate 
these salient features of financial returns distributions, namely more probability mass in 
the tails and around the center of the distribution than in the benchmark Gaussian case. 
Capturing these features of financial markets is crucial for many purposes including but 
not limited to: derivatives pricing, hedging, forecasting volatility, portfolio management, 
regulatory issues, value- at- risk, and so on. 

The traditional ARCH specification expresses the current volatility level as a function 
of past squared shocks (about a mean or average return) at various lags. This implies 
that the volatility is not constant over time anymore but that it depends on how large 
recent deviations (positive or negative) from the mean have been. In the more general 
GARCH specification, the current volatility level is still a function of past squared shocks 
to the returns, but it is also a function of past lagged levels of itself, thus making the 
model even more flexible. The GARCH(p,q) specification for the volatility of at time t 
can be described by the following equation: 



af =a a +J2 OiA-i + U Pi°t-i (2°) 

i=l z=l 

where q and p are the maximum lags allowed by the model for past shocks and volatility 
levels respectively, and tt-i is the return shock (about the mean) at lag i. 

One can obviously allow q and p to be as large as one wants, but a GARCH(1,1) 
model in practice turns out to be surprisingly flexible. Moreover, overfitting is often 
the recipe for poor out-of-sample performance, and parsimonious models often end up 
defeating more complex ones when tested outside of the in-sample period. Hansen and 
Lunde [26| compare 330 (G)ARCH-type models in terms of their ability to describe the 
conditional volatility of exchange rate and IBM return data, and find no evidence that 
a GARCH(1,1) is outperformed by more sophisticated models in the analysis. Therefore 
we choose to focus on the parsimonious, yet very apt, GARCH(1,1) model for purposes 
of comparing the dynamics of time- varying volatility in our lattice-generated returns and 
in NASDAQ historical returns. Our model can thus be written as 

o\ = a + ai^_x + h^t-\ . (21) 

The goal here is twofold. First, it is to investigate whether our lattice model is able 
to produce returns volatility dynamics displaying some form, if any, of ARCH/GARCH 
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effects. Second, it is to examine whether our lattice model is able to produce returns 
volatility dynamics that are rather consistent with those of NASDAQ historical returns. 
In figures ITU1 through IT51 for each set L and N (lattice and NASDAQ returns), we fit a 
GARCH(1,1) model onto the returns data through a Levenberg-Marquardt optimization 



Before getting to the tables, one may first notice that figures ITU1 through [Hi! display 
similar-looking sets of charts when one compares the lattice-generated graphs with the 
NASDAQ historical ones. The top portion represents returns (gains or losses) over time, 
and one can readily see that there are clusters of volatility in each, indicating that the 
model seems to be capable of capturing them. The middle portion simply represents 
the evolution of the asset value over time. Visually, the price dynamics generated by 
the lattice model appear 'plausible' as those of a financial market such as the NASDAQ. 
Finally, the bottom plot represents the evolution of the volatility over time, for both 
the lattice and the NASDAQ data. Here again, the volatility patterns generated by the 
lattice model seem credible as those of financial markets, displaying a variety of spikes 
of various sizes and frequencies. 

Underneath the plots, each table displays the estimated parameters ao, a\ and 
followed by their standard errors in parenthesis and their t-statistics, in square brack- 
ets, for both the lattice and the NASDAQ returns data. Except for the NASDAQ pa- 
rameter ao in Tab. 2J every parameter is statistically significant at the 5% significance 
level, and most parameters are even statistically significant at the 1% significance level. 
This indicates that both the NASDAQ returns and our lattice-generated returns display 
ARCH/GARCH behavior in volatility. Moreover, and more remarkably, the estimated 
parameters coming from our lattice-generated returns are extremely close, especially in 
magnitude, to the estimated parameters coming from the NASDAQ historical returns. 
For instance, in Tab. [H ao, a\ and /?i are estimated to be 4.14 x 10~ 9 , 0.022450 and 
0.966177 respectively for the lattice-generated model, while they are 4.20 x 10 -9 , 0.031982 
and 0.951740 respectively for the NASDAQ historical returns. Although this does not 
indicate forecasting abilities on the part of the lattice model, it does show that it is able 
to reproduce, with precision, important financial markets volatility dynamics and thus 
has the potential to provide future insights on the inner mechanics of such markets. 

5. Summary and conclusion 

By their very nature, historical market data constitute only instances of some stochas- 
tic process. It is thus desirable to have available a stochastic model of a financial market. 
As a result, such a model grants complete access to the market's stochastic features 
through the act of drawing multiple instances, or ensembles. This allows the possibility 
of a detailed investigation of market dynamics and the features that define them. 

We have studied the properties of a stochastic lattice model which by design derives 
its features from next-nearest-neighbor interactions of microscopic entities that live on 
a linear chain in discrete time. The model is realized as a lattice field theory with field 
components being interpreted as market returns, and is subject to computer simulation 
with an updating algorithm inspired by the evolution model by Bak, Tang and Wiesenfeld 
Q that drives the lattice field into a self-organized critical state. We present evidence, 
including power law behavior, that the critical state is indeed achieved. Its presence is 
the salient feature of the model. 
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We compute time series of market returns, prices, volatilities, and returns frequency 
distributions, all of which are remarkably consistent with historical market data for the 
NASDAQ stock index. In particular the 'fat tails' feature of the returns distribution 
comes out effortlessly. It is worth noting that, aside from (trivial) units, the lattice model 
has no adjustable parameters. Perhaps most remarkable is the observation that standard 
financial industry analysis tools, in our case the GARCH(1,1) model [IBj], produce fits that 
make the lattice model time series almost indistinguishable from real financial market 
data. 

It seems that the most important conclusion derived from our study is that self- 
organized criticality, being a fundamental driving force for the model, is also key to 
characterizing real world financial markets. To say the least, self-organized criticality 
should be an important component of intrinsic market dynamics, allowing us to model 
financial instruments using lattice methods which may hopefully be competitive with 
current industry practice. 

Of course, self-organization is surely not the only mechanism that drives financial 
market dynamics. Among other things, trading occurs within a background of arbitrage 
opportunities [3[ • We are looking forward to combining past studies focused on arbitrage 
[2| with the present model in order to develop an even more realistic stochastic market 
model. 
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